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Abstract 

We develop an idea of Evans and O'Connell [11], Englander and Pinsky [8] and 
Duquesne and Winkel [3] by giving a pathwise construction of the so called 
'backbone' decomposition for supercritical superprocesses. Our results also com- 
plement a related result for critical (1 + /3)-superprocesses given in Etheridge 
and Williams [9]. Our approach relies heavily on the use of Dynkin-Kuznetsov 
N-measures. 
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1 Introduction 

In Evans and O'Connell [11], and later in Englander and Pinsky [8], a new decom- 
position of a supercritical superprocess with quadratic branching mechanism was in- 
troduced in which one may write the distribution of the superprocess at time t > 
as the result of summing two independent processes together. The first is a copy of 
the original process conditioned on extinction. The second process is understood as 
the aggregate accumulation of mass from independent copies of the original process 
conditioned on extinction which have immigrated 'continuously' along the path of an 
auxilliary dyadic branching particle diffusion which starts with a Poisson number of 
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particles. The embedded branching particle system is known as the backbone (as op- 
posed to the spine or immortal particle which appears in another related decomposition, 
introduced in Roelly-Coppeletta and Rouault (1989) and Evans (1993)). In both [TT] 
and [8] the decomposition is seen through the semi-group evolution equations which 
drive the process semi-group. However no pathwise construction is offered. 

In Duquesne and Winkel [3] a version of this decomposition which, albeit does not 
take account of spatial motion, was established in much greater generality. In their 
case, quadratic branching is replaced by a general branching mechanism ip which is the 
Laplace exponent of a spectrally positive Levy process and which satisfies the condi- 
tions < —ip'(0+) < oo and f°° l/^(£)d£ < oo. Moreover, the decomposition is offered 
in the pathwise sense and described through the growth of genealogical trees embedded 
within the underling continuous state branching process. In their case the backbone 
is a continuous-time Galton Watson process and the general nature of the branching 
mechanism induces three different kinds of immigration. Firstly there is continuous 
immigration which is described by a Poisson point process of independent processes 
along the backbone where the rate of immigration is given by a so-called excursion 
measure which assigns zero initial mass, and finite life length of the immigrating pro- 
cesses. A second Poisson point process along the backbone describes the immigration of 
independent processes where the rate of immigration is given by the law of the original 
process conditioned on extinction and with an initial mass randomised by an infinite 
measure. This accounts for so-called discontinuous immigration. Finally, at the times 
of branching of the backbone, independent copies of the original process conditioned 
on extinction are immigrated with randomly distributed initial mass which depends on 
the number of offspring at the branch point. The last two forms of immigration do not 
occur when the branching mechanism is purely quadratic. 

Concurrently to the work of j3] and within the class of branching mechanisms corre- 
sponding to spectrally positive Levy processes with paths of unbounded variation (also 
allowing for the case that —ip'(0+) = oo), Bertoin et al. [2] identify the aforementioned 
backbone as characterising prolific individuals within the genealogy of the underling 
continuous state branching process. Here, a prolific individual is understood to be an 
individual whose descendants become infinite in number. 

In this paper we develop the decomposition of Duquesne and Winkel [3] further by 
adding in the following additional features. We allow the possibility J 1/V>(£)d£ = oo 
which includes the possibility of supercritical processes whose total mass may, with 
positive probability, die out without this ever happening in a finite time. This also 
allows the inclusion of branching mechanisms which belong to spectrally positive Levy 
processes of bounded variation (previously excluded in [2] and [3]). Secondly our de- 
composition takes care of spatial motion of individuals, thereby bringing the Duquesne- 
Winkel decomposition back into the setting of sup erpro cesses. Finally, in the case that 
we ignore spatial motion, our analysis also allows for the case that — ^'(0+) = oo. Our 
proof is fundamentally different to that of [3] and relies largely on the manipulation 
of the semi-group evolution equations in the spirit of [TT], taking advantage of the so 
called N- measure of Dynkin and Kuznetsov [7]. 

The remainder of the paper is structured as follows. In the next section we introduce 
some preliminary notation and remind the reader of some standard results relevant to 
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the subsequent exposition. In Section [3~T1 we describe a branching particle diffusion on 
which independent superprocesses immigrate in three different ways. In particular we 
give a key result in which the semi-group of the aforementioned process with immigra- 
tion is characterised. With the latter in hand, we are able to state and prove in Section 
I3.2[ the backbone decomposition for supercritical superprocesses. Finally in Section H] 
we give a proof of the the key analytical result in Section 13.11 Along the way we shall 
also establish the slightly stronger backbone decomposition for the case of continuous 
state branching processes (ie. when spatial considerations are ignored). 



2 Preliminaries 

In this section we outline some standard notation and mathematical tools as well as 
key existing results, all of which will be the ingredients that together will make up the 
main result. 



2.1 (V, ^)-superprocess 

Suppose that X = {X t : t > 0} is any superprocess motion on R d which is well defined 
for initial configurations in .MfQ^), the space of finite and compactly supported mea- 
sures, having an associated conservative diffusion semi-group V := {Vt ■ t > 0} on R d 
and general branching mechanism ip taking the form 

ip(\) = a\ + (3\ 2 + [ (e~ Xx - 1 + \xl {x<1} )n{dx) 

for A > where a G R, > and II is a measure concentrated on (0, oo) which 
satisfies f, Q ^(1 A x 2 )U(dx) < oo. This implies that the total mass of the process X is 
a continuous-state branching process with branching mechanism ijj for which standard 
references, e.g. [5j [10], dictate that we need to assume that —ip'(0+) < oo. Note 
however that without this condition it is always the case that — ^'(0+) G (— oo, oo] and 
within this regime, continuous-state branching processes are always well defined; see 
for example [2]. To exclude the case of explosive behaviour, we assume throughout 
that 




d£ = oo. 



Moreover, we insist that if) (oo) = oo which means that with positive probability the 
event lim^oo | \X t \ \ = will occur; see for example the summary in Chapter 10 of 
Kyprianou [T7]. We refer to such processes throughout as (V, -^-superprocesses. 

Remark 1. It is worthy of note that the assumption that V is a conservative diffusion 
semi-group on R d can easily be replaced throughout by the much weaker assumption 
that V is a general Borel right Markov process with Lusin state space, just as in [11] 
or [U Ej, at no cost to the analysis. Indeed all of the proofs go through verbatim. 
However, purely for the sake of presentation, we keep to the more familiar Euclidian 
setting. 
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Remark 2. Whilst the vast majority of all literature concerning (V, ^)-superprocesses 
requires the branching mechanism satisfies — ip'(0+) < oo, an example for which an 
infinite branching rate is permitted can be found in Fleischmann and Sturm [12] . 

For each fi G A4p(M. d ), we denote the law of X with initial configuration /x by P M . 
The following standard result from the theory of superprocesses (c.f. Dynkin [U [5] for 
example) describes the evolution of X as a Markov process. 

Lemma 1. For all f G bp(R d ), the space of non-negative, bounded and measurable 
functions on M. d , 

-logE M (e~ (/ ' Xt> ) = / u f (x,t)fi(dx), fi G M F (R d ), t > 

where Uf(x,t) is the unique non-negative solution to the integral equation 

u f (x,t)=V t [f](x)- f\-V,[rl>{u f {;t-s))}{x). (2.1) 

Jo 

Here we have used the standard inner product notation, for / G bp(M. d ) and \i G 
J\4(M. d ), the space of measures on M. d . 

(/,/*>= / /(x)/i(dx). 

Accordingly we shall write ||/i|| = 

Remark 3. In the case that we take P to correspond to a particle remaining stationary 
at a point, equation (12. ip collapses to the classical integral equation describing the 
evolution of a continuous state branching process. As alluded to above, it is known in 
this case that a unique non-negative solution exists, even in the case that — ip'(0+) = 
— oo. 

2.2 Criticality 

As noted above the total mass of a (V, ?/;)-superprocess is a continuous-state branching 
process with branching mechanism ijj. Since there is no interaction between spatial 
motion and branching we can therefore characterise the (V, "0)-superprocess into the 
catagories of supercritical, critical and subcritical accordingly with the same catagories 
for continuous-state branching processes. Respectively, these cases correspond to ?//(0+) < 
0, ip'(0+) = and ^'(0+) > 0. Recall that even when X is supercritical, it is possible 
that the process becomes extinguished, i.e. liim.-^ \ \X t \\ = 0. The probability of the 
latter event is described in terms of the largest root, say A*, of the equation = 0. 
Note that it is known (cf. [H]) that ip is strictly convex with ip(0) = and hence since 
if) (oo) = oo and ^'(0+) < it follows that there are exactly two roots in [0, oo), one of 
which is A* and the other is 0. For \x G A4p(M. d ) we have 

P M (lim||Xt|| = 0) = e~ x * M . (2.2) 

tfoo 
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We also recall that, if in addition 

then the event {lim^oo \ \X t \ \ = 0} agrees with the event of extinction, namely {( < 00} 
where 

C = inf{t > : ||X t || = 0}. 

Moreover, when the integral test in (12. 3p fails, the supercritical continuous-state branch- 
ing process becomes extinct with zero probability. This means that the event of be- 
coming extinguished corresponds to the total mass trickling away to zero but none the 
less being strictly positive at all finite times. An example of an infinite mean supercrit- 
ical branching mechanism for which the phenomena of becoming extinguished but not 
extinct is = X — \ a where a G (0, 1). A second example in this class is Neveu's 
branching mechansim ^(A) = A log A. Note that in the first example, the associated 
spectrally positive Levy process has paths of bounded variation as it can be written 
as the difference of an a-stable subordiantor and a linear unit-rate drift. The second 
example corresponds to a spectrally positive Levy process with paths of unbounded 
variation as in that case it is known that the underlying Levy measure is given by 
n(dx) = x~ 2 dx. 



For the remainder of the paper, unless otherwise stated, we shall 
henceforth assume only that 'ip is a non-exploding, supercritical branching 

mechanism satisfying —'ip'(0+) < 00. 

It is well known that there is a link between ip and another branching mechanism 
ip* where, for A > —A*, 

:= ^(A + A*) 

= a*X + (3\ 2 + [ (e~ Xx - 1 + Xxl {x<1} )e- x * x U(dx) (2.4) 
J (0,00) 

and 

a* = a + 2/3\*+[ (1 - e~ x * x )xll(dx). 

J (0,1) 

The connection between ip and i/j* has a distinct probabilistic interpretation that we 
shall now briefly discuss. 

Recall that a continuous-state branching process with branching mechanism ip can 
always be written as a time-changed Levy process with no negative jumps and whose 
Laplace exponent is precisely ijj; see for example [18j [17]. The assumption ip'(0+) < 
implies that the underlying Levy process drifts to +00. Moreover, by a classical result, 
the branching mechanism ip* corresponds to the underlying Levy process conditioned 
to drift to —00 (see Exercise 8.1 of [TT] ) . As the next lemma confirms, it also turns 
out that ift* is the branching mechanism of a superprocess which can be identified as 
the (V, ip) -superprocess conditioned to become extinguished. (Similar results can be 
found in [2], Abraham and Delmas |T| and Sheu [21]) 
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Lemma 2. For each \x G A4_p(M ), define the law of X with initial configuration fi 
conditioned on becoming extinguished byF*^ (with expectation operator K* ). Specifically, 
for all events A, measurable in the natural sigma algebra of X , 

P;(A)=P,,(A| Jim 1 1^11 = 0). 

Then, for all uniformly bounded f : M. d — >■ [—A*, oo), 

-logE;(e~ (/ ' Xi) ) = f u* f (x,t)/i(dx) 

where 

u* f {x,t)=u f+x *(x,t)-\* (2.5) 
and it is the unique solution of 

u}(x,t) = V t [f](x)- [ ds-P s [r(u* f (;t-s))](x) (2.6) 

Jo 

where = + for A > —A*. That is to say (X, P*) is a (V , tp*)-superprocess. 

Proof. Define the event £ = {lim^o \ \Xt\\ = 0}. Making use of the Strong Markov 
property and (12. 2\\ , we have for / G bp(R d ), 

E *( e -</,^>) = E,(e-^\S) 

= e x *^%(e-^e- x *^) 

It is trivial to check that w/+a* (•, t) — A* solves (I2.6p . Moreover, since ip* is the Laplace 
exponent of a spectrally positive Levy process and ^*'(0+) = ^>'(A*) > 0, it follows 
that the solution to (12. 6 p is unique by Lemma [TJ □ 

2.3 N*-measure 

Associated to the laws {P| : x G M. d } are the measures {N* : x G defined on the 
same measurable space, which satisfy 

K(l - e-CW) = -logE; a (e-^>) (2.7) 

for all / G bp(M. d ) and t > 0. Such measures are formally defined and explored in 
detail in [7]. The measures {N* : x G M d } will play a crucial role in the forthcoming 
analysis. Intuitively speaking, the branching property implies that P| is an infinitely 
divisible measure on the path space of X, X := JH(M. d ) x [0, oo), and (12. 7p is a 'Levy- 
Khinchine' formula in which N* plays the role of its 'Levy measure'. In the context of 
[3], the measure N* is the analogue of what Duquesne and Winkel call the excursion 
measure, however, whilst the latter encodes genealogical trees, N* does not. 
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2.4 Prolific individuals 

In Duquesne and Winkel [3] and Bertoin et al. [2] it was shown that there are certain 
geneaologies embedded in supercritical continous state branching process which are 
exclusively responsible for the infinite growth of the process. They show that one may 
identify such geneaologies in the form of a continuous-time Galton Watson process (that 
is to say, a version of the Galton Watson process in which individuals remain alive for 
an independent and exponentially distributed period of time with a common rate before 
splitting). The generator of such a continuous-time Galton Watson processes is usually 
identified in the form 

F(s) = qJ2Pn(s n -s), 

n>0 

where q > is the common rate of splitting and {p n : n > 0} is the offspring dis- 
tribution. For the particular continuous-time Galton Watson process representing the 
prolific geneaology in a supercritical continuous-state branching process with branching 
mechanism ip, the aforementioned authors show that 

F(s) = ^(\*(l-s)),se(0,l). (2.8) 

Moreover the individual components of F are given by q = po = pi = and for 

n > 2, 

p " = {^*) 2l <-> + (*rJ M ^ n(dx)} . (2.9) 

Duquesne and Winkel [3] go further and show that, when ip'(0+) £ (— oo,0) and 
(12. 3p holds, the law of a continuous state branching processes with branching mecha- 
nism i[) is equal to that of a process in which immigration occurs on the continuous-time 
Galton Watson process of prolific individuals in three different ways. These are, two 
types of Poisson immigration along the life span of each prolific individual and an 
additional package of immigration at each point of fission of prolific individuals. In 
the latter case, if a prolific individual has n prolific offspring then a continuous-state 
branching process with branching mechanism ip* immigrates at that moment of time 
with random initial mass given by the distribution 

Vn(dy) = M ^, (A ^ {/3(A*) 2 5 (dy)l {n=2} + (ATge-^n(dy)} ■ (2.10) 

In the next section we progress the result of [3] by relaxing their assumptions on if; 
to include the cases that ^'(0+) = — oo and f°° 1/V(£)d£ = oo as well as taking 
into account spatial considerations at the expense of keeping the condition —ip'(0+) £ 
(0,oo). 
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3 Backbone decomposition 



3.1 A branching particle diffusion with three types of immi- 
gration 

Let A^ a (M d ) C A^i?(R d ) be the space of finite atomic measures on M. d . We shall write 
Z for a branching V- motion whose total mass has generator given by (12. 8p . Hence Z 
is the M. a (M, d )- valued process in which individuals, from the moment of birth, live for 
an independent and exponentially distributed period of time with parameter 
during which they execute a "P-diffusion issued from their position of birth and at 
death they give birth at the same position to an independent number of offspring with 
distribution {p n : n > 2}. We shall also refer to Z as the (V, F)-backbone. Its initial 
configuration is denoted by v £ Ai a (M. d ). Moreover, when referring to individuals in 
Z we may use of classical Ulam-Harris notation, see for example p290 of Harris and 
Hardy [15]. The only feature we really need of the the Ulam-Harris notation is that 
individuals are uniquely identifiable amongst T, the set labels of individuals realised 
in Z. For each individual u £ T we shall write t u and a u for its birth and death times 
respectively and {z u {r) : r £ [t u , cr u ]} for its spatial trajectory. 

Inspired by [3] and [2], we are interested in immigrating independent (V,^*)- 
superprocesses on Z in a way that the immigration rate is related to the subordinator 
whose Laplace exponent is given by 

0(A) = ^*'(A) - = + A*) - (3.11) 

together with some additional immigration at the splitting times of Z. Note in partic- 
ular that the right hand side of (13. lip can be written more explicitly in the form 

0(A) = 2/3A + f (1 - e- Xx )xe- x * x Il(dx). 

</(0,oo) 

Definition 1. For v £ M..JW 1 ) and /i £ M. F {R d ) let Z be a (V, F) -branching diffusion 
with initial configuration v and X an independent copy of X under P* Then we define 
the measure-valued stochastic process A = {A t : t > 0} on~R d by 

A = X + I w + I ¥ * + P 

where the processes I N * = {if : t > 0}, I F * = {If : t > 0} and P> = {I? : t > 0} are 
independent of X and, conditionally on Z , are independent of one another. Morover, 
these three processes are described pathwise as follows. 

(i) Continuous immigration: The process I N * is measure-valued on M. d such that 

ueT t At u <r<tAcr u 

where, given Z , independently for each u £ T such that t u < t, the processes 
^(i,«,0 are coun f a })i e { n number and correspond to X-valued, Poissonion immi- 
gration along the space-time trajectory {(z u (r),r) : r £ (r u ,t A a u ]} with rate 
2/?dr x dW Zu{r) . 
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(ii) Discontinuous immigration: The process I v * is measure-valued on M d such that 

IT := £ £ 

U&T tf\Tu<r<t/\<J n 

where, given Z , independently for each u G T swc/i t/iai r u < £, i/ie processes 
j^(2,«,r) countable in number and correspond to X-valued, Poissonion immi- 
gration along the space-time trajectory {(z u (r),r) : r G (T u ,tAer u ]} with rate 

(iii) Branch point biased immigration: The process I v is also measure valued on 
R d such that 

TV -1 r , y(3,M) 

where, given Z, independently for each u G T siic/i £/iat a u < t, £/ie process X^' u ^ 
is an independent copy of X issued at time a u with law Fy u s x , where Y u is an 
independent random variable with distribution r)N u (dy). 

Moreover, we denote the law of A by P^. 

Remark 4. In the very special case that ip(\) = — aA + bA 2 , where a, b > 0, note 
that the discontinuous and branch point biased immigration are absent. Moreover, 
= aA + bA 2 , the backbone has binary splitting and therefore agrees with the 
backbone in Evans and O'Connell jllj . 

Note that the total mass Z t (M. d ) of the backbone is the continuous-time Galton- 
Watson process of prolific individuals found in Bertoin et al. j2]. Note also that the 
process ((A, Z), P^xu) is Markovian. This is immediate from three important facts. 
Firstly the backbone, Z, is a Markov branching diffusion. Secondly, conditional on 
Z immigrating mass occurs independently according to a Poisson point process or as 
additional indpendent packages at the splitting times of Z. Finally, the mass which 
has immigrated by a fixed time evolves in Markovian way thanks to the branching 
property. Indeed, using these facts it is not difficult to justify that, for all s, t > 

E MXi/ (e-^ A *+^|{(A u ,^) :u<t}) = h.(A t ,Z t ), 

where for m G M{R d ), n G M a {R d ) and s > 0, h s (m,n) = E mxn (e" (/ ' As) ). 

We conclude with the main result of this section which, amongst other things, shows 
that A is a conservative process. The proof is given in section HI 

Theorem 1. For every /i G M. F (R d ) } v G M a (M d ) and f,h G bp(M. d ) we have 

E^ xu (e- {f ' At) - {h ' Zt) ) = e -W'> t )>i*)-< v fM-fi>'') (3.12) 

where exp{— Vf^(x, t)} is the unique [0, l]-valued solution to the integral equation 

e -v f M*,t) = 7? t [ e -^]( x ) + _L /" < ds .p s [^*(_A*e-^ h( -' t - s) +n}(-,t- S ))-^(nK-,t- S ))](x). 

A * Jo 

(3.13) 

for x G M. d and t > 0. In particular, for each t > 0, A t has almost surely finite mass. 
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3.2 Prolific backbone decomposition of a supercritical (V, ifj)- 
superprocess. 

A consequence of Theorem [1] is the following theorem which constitutes our main result. 
It deals with the case that we randomise the law P MX f for fi G A4,F(IR d ) by replacing 
the deterministic choice of v with a Poisson random measure having intensity measure 
X*fi. We denote the resulting law by P M . 

Theorem 2. For any \i G M.F(M. d ), the process (A, P M ) is Markovian and has the 
same law as (X, P M ). 

Proof. The proof is guided by the calculations found in the proof of Theorem 3.2 of [UJ . 
We start by addressing the claim that (A, P M ) is a Markov process. Given the Markov 
property of the pair (A, Z), it suffices to show that given A f , the atomic measure Z t is 
equal in law to a Poisson random measure with intensity X*A t . Thanks to Campbell's 
formula for Poisson random measures (see e.g. Section 3.2 of [IE]), this is equivalent 
to showing that for all h G bp(JBL d ), 

E,(e-^\A t ) = exp{-(A*(l - e" h ), A t }}, 

which in turn is equivalent to showing that for all f,h E bp(JBL d ), 

E/i ( e -</,A.Mft,s*>) = E^e-W-'^+'M). (3.14) 

Note from (13.121) however that when we randomise v so that it has the law of a Poisson 
random measure with intensity A*/x, we find the identity 

E ( e -(fM)-(h,Zt)j = e -(u*(;t)+\*(i-e- v fM-^)^)^ 

Moreover, if we replace / by A*(l — e~ h ) + / and h by in (I3.12p and again randomise 
v so that it has the law of a Poisson random measure with intensity A*/i then we get 

E/i ( e -<A*(l-e- ft )+/,At>) = e -K, (1 _ e - h)+/ (-. i )+ A *( 1 - cx P{-^.(i- e -'>) + /,o(-' i )})^) _ 

These last two observations indicate that (I3.14p is equivalent to showing that for all 
f,he bp(R d ), iGl d and t > 0, 

u}(x, t) + A*(l - e~ v ^ l) ) = u^ {1 _ e . h)+f {x, t) + A*(l - e -^(i-e^)+/,o(^)). (3.15) 

Note that both left and right hand side of the equality above are necessarily non- 
negative given they are the Laplace exponents of the left and right hand sides of (13.141) . 
Making use of (12.61) and (13.131) . it is computationally very straightforward to show that 
both left and right hand side of (J3.15p solve (12. ip with initial condition / + A*(l — e~ h ). 
Since ( 12. ip has a unique solution with this initial condition, namely Uf + x*n- e -h\(x,t), 
we conclude that ( 13 . 1 5 j) holds true. The proof of the claimed Markov property is thus 
complete. 
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Having now established the Markov property, the proof is complete as soon as we 
can show that (A, P^) has the same semi-group as (X, P M ). However, from the previous 
part of the proof we have already established that when f,hE bp(R d ), 

E„(e~ {f ' At) ~ {h ' Zt) ) = e - <u A*(i- e -")+/^) = E (e~ {f+x * {1 ~ e ~ h) ' Xt} ). 

In particular, choosing h = we find 

E,(e-^ A <>)=E,(e-^) 

which is equivalent to the equality of the semi-groups of (A, P M ) and (X, P M ). □ 

Remark 5. In the proof above, we have established the so-called Poissonisation prop- 
erty of superprocesses and continuous state branching processes. Namely that, when 
treating \*X t as an intensity measure of a Poisson random field, one generates a set of 
points whose positions are equal in law to the support of Z t . Fleischmann and Swart 
[13] appeal directly to this idea to analyse the law of Z in terms of X. 

Remark 6. Once the reader is familiar with the main ideas of Theorem [2] it should 
be quite clear how to describe in a pathwise sense the backbone-type decomposition 
in Englander and Pinsky [8]. In their paper, they work with a spatially dependent 
branching mechanism A) = a(-)A + b(-)A 2 . Given the semi-group computations in 
[8] one may easily construct the associated pathwise decomposition. There is no dis- 
continuous immigration and no branch point biased immigration. However continuous 
immigration does occur along the backbone at rate 2b(-)d£ x dN* where N* is again 
the measure constructed in [7] which is related to law of the superprocess conditioned 
on extinction. The latter, as well as the law of the backbone, are already described in 
analytical detail in [8]. 

3.3 Prolific backbone decomposition of a supercritical continuous- 
state branching process 

The analysis leading to the proof of Theorem [1] also reveals that the assumption that 
—ip'(0+) < oo can be dropped when considering the backbone decomposition for con- 
tinuous state branching processes. Formally we state this as a theorem. 

Theorem 3. When V corresponds to a particle remaining stationary at a point, say 
0, the conclusion of Theorem [H still holds for all fi = xSq with x > 0, even when 
-7/7(0+) = oo. 



4 Proof of Theorems [T] and [3 

To prove Theorem [1] it suffices to show, thanks to Lemma [2j that for all f,hE bp(M. d ), 
v e .M a (M d ) and t > 0, 

E^^e-tf' 1 **-^) = e -WA(-A"> (4.I6) 

where / := I N * + I F * + I v and Vf t h solves (I3.13p . We do this with the help of some 
preliminary results. 
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Lemma 3. For all f 6 bp(R d ), v e M a (M. d ), fi e M F {R d ) and t > 0, we have 

E, xu (e-^* +I t^\{Z s :s<t}) = exp [~ j\<t> ° «/(', * - *), Z s >ds} 
Proof. Using the notation from Definition [U write 

</,/r = £ £ a^ r) )+E E </.*£r u,r) >- 

u£T tAr u <r<tAcr u u£T tAr u <r<tAcr u 

Hence conditioning on Z, appealing to independence of the immigrating processes 
together with Cambell's formula and (j2.7p we have 

E^ u (e-^* +I t^\{Z s :s<t}) 

{rtAvu 
-£2/3 / dr-N^ (r) (l-e-^->) 

rtAau r 

-E / / dr x ^~ A * yn ( d ^) • E k„ w c 1 - e ~ (/,Xt ~ r: 

ue -r^ tATu ^ (°'°°) 

{rtAa u 
~£ 2 ^ / dr ' u *f( z u( r ),t - r) 

rtAcF u /> 

-E / / drxj/e-^nfd^.a-exp^^^.t-rM) 
ue Y tATu ^(°>°°) 

= expj-^ dr- (0o M *(-,t-r),Z r )| 

as required. □ 

Lemma 4. Suppose that f,h G op(IR d ) and <7 s (x) jointly measurable in (s,x) and 
bounded on finite time horizons of s. Then for x G lR d and £ > 0, 

E UJt „ fexp |- jf (o t _ s , Z s )d s - (/, i?> - (/i, Zt) Y) = e-WM 



where exp{— u(x, £)} i/ie unique [0, l]-valued solution to the integral equation 
e -^,t) =Vt[e -h ]{x) + ^ jy s . Vs[Ht ^ (4 _ 17) 

and, /or A > —A*, 

Ht-,(;\) := A^'(A*) + (3\ 2 + / (e~ A:l -l + Ax)e~ {A * +u / ( ' i - s))a; n(dx). 
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Proof. Following similar arguments to those in the proof of Theorem 2.2 of Evans and 
O'Connell [UJ, it suffices to consider the case that, in addition to the assumptions 
in the statement of the lemma, g is time- invariant. Moreover, using the branching 
property of Z it suffices to consider the case that v = 5 X for x £ M. d . In that case, 
suppose that £ := {£ t : t > 0} is the stochastic process whose semi-group is given by V . 
We shall use the expectation operators {E x : x £ W 1 } defined by E x (f(£ t )) = V t [f](x). 
Define a new semi-group (of the diffusion £ killed at rate g) 

V?[f](x) = E x (e-ti^ ds f(£ 



for f,g £ bp(M. d ). Standard Feynman-Kac manipulations (see Lemma 2.3 of [TT]) give 
us that 

=*>*[/](*)- fds-VM-^UfMix). (4.18) 
Jo 

Conditioning on the time of the first branching and recalling that branching occurs at 
rate q = ip'(X*) we get that 



j(x,t) 



= e~ qt Vf[e- h ](x) + q [ ds ■ e~ qs V 9 s J^P^'^''^ I Vn(dy)e- yu ^' t - s) 

J0 ' [ n > 2 7(0,00) 

Next note from ( 12~TU|) that 
Y.Pne-™^-* [ V n(dx)e- xu f^- s) 

n>2 J(0,oo) 

V/3(A*e- w( -' i " s) ) 2 l {n=2} + — [ (xX*e- UJ ^ t - s) ) n e~ {X ' +u 'f { '' t ~ s))x Il(dx) 

n>2 U - J ^oo) 



(x). 



1 

gA* 
1 

q\* 



f3(X*e 



-«(-,t-a)\2 



'(0,oo) 



- 1 - x\*e- uj{ -> t - s) )e- {x * +u *f { -> t - s))x Tl(dx) 



(0,oo) 



-^-[#t-.(-, -A*e-" ( -'*- s) ) + gA'e - "*'*-)] 



r/A 

We now have that 

-u(x,t) 



1 /** 

= e _9t ^f[e- ft ](z) + — / ds • e^Pf -A^e^^'"^) + gA*e^ ( -'*- s) ] (x) 

A* Jo 

= V?[e~ h ](x) + ^£ds- VI [H t _ s (; -A*e^<^)] (x) (4.19) 

where the second equality follows by a standard technique found, for example, in 
Lemma 4.1.1. of [5]; see also the computations in [TT] . 
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Next, we use ( I4.18P and note that 



-uj(x,t) 



V t [e- h \{x)- f ds-V s [g(-)Vl s [e- h ](-)}(x) 
Jo 

+1 jf ds ■ {v s [H t „ s (; -A*e-"(-'*-))](x) 

- ^ dr • P r ^(.)Pf_ r [£Tt- s (-, -A* e -^^))]](x)} 



A* . 

where in the final equality we have used (I4.19p to deduce that 

I'd, • V s [g{-)Vl s [e- h \{-)\{x) + ±-j\ s -j\r- V r [g(-)V 9 s _ r [H t -s(; -A*e"^^)]](x) 



* 1 



ds-V s [g(-)VUe- h ](-)](x) 



X* 



dr-V r 



&s-V s [g{-)Vl s [e- h ]{-)]{x) + - \ dr-V r 



' dd- V 9 e [H t _e- r (; -AV-*< '•*-^)](-) 



(*) 



dr ■ P r 



(•) + 



1 



i-r 



dd ■V 9 e [H M (.,-\*e-^> t - r -V)}(-) 



(x) 



f ds-V s [g( 
Jo 



-u}(-,t—s) 



X . 



The proof is complete as soon as we can establish uniqueness to (14.1 7ft . By mul- 
tiplying the latter equation through by A* we note that, by an application of Lemma 
2.1 of [11] (which offers sufficient conditions for solutions to a general family of inte- 
gral equations), it has a unique solution providing the assumptions of that lemma are 
satisfied. For this purpose it suffices to check that for each y £ M. d and A £ [0,A*], 
J(s, y, A) := [H s (y, —A) — g(y)X\ is continuous in s and that for each fixed T > 0, there 
exists a K > such that 

sup sup | J(s,y,u(y)) - J(s,y,v(y))\ < K sup \u(y) - v(y)\ 

where u and v are any two measurable mappings from R d to [0, A*]. In light of the 
assumption of boundedness on g(y), thanks to the triangle inequality, it suffices to 
check that for each fixed T > 0, there exists a K > such that 

sup sup \H 8 (y, -u(y)) - H s (y, -v(y))\ < K sup \u(y) - v 

where u and v are any two measurable mappings from R d to [0, A*]. 
To this end let us define for A > —A* and u > 0, 



Xu{\) := A^'(A*) + /3A 2 + 



f (e~ Xx - 1 + \x)e~ {x * +u)x Il(dx). 



r (0,oo) 
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so that by definition H s (y, A) = Xu*( y ,s)W, for A > —A*. We need the following facts 
about XuW'- 

Lemma 5. For A > —A* we have that 

xu(a) = v*(a + u)- r{u) - \\r\u) - r'm- 

Moreover, when we allow ^'(0+) G [—00, 00), for each u,u>0 we have 

sup_ sup |xl(-A)| < 00. 

u<u<u AG[0,A*] 

If however, —ip'(0+) < 00 then we may take u — in the above inequality. 

The proof of this result is somewhat technical and disjoint from the core of the 
argument that we are currently pursuing, so its proof is postponed until the end of the 
paper. 

With the help of the above lemma, we see that for each fixed T > 0, 
sup sup \H s (y, -u(y)) - H s (y, -v(y))\ 

= SUp SUp \Xu*Jy,s)(-u(y)) ~ Xu* f (y,s)(-V(y))\ 

< sup_ sup \Xu*(-u(y)) - Xu*(-v(y))\ (4.20) 

0<m*<m t y£R d 

< K sup \u{y) - v{y)\ 

y£R d 

where u and v are any two measurable mappings from R d to [0, A*], 

K = sup sup |Xu*(— A)| < 00 (4-21) 

0<tt*<u T Ae[0,A*] 

(observe that the inequality is true if and only if ^'(0+) > —00) and 

ut = sup sup u*f(y, s) < 00. 

Note that the finiteness of ut can be deduced as follows. Suppose, without loss of 
generality, that / is bounded by 9 > 0. Then for all y G M. d and s > 0, 

e -u}(y,s) = E y e -</,**>) > E^(e~ e||Xs|1 ) = e~ u ^ s) 
where Ug(s) is the unique solution to the equation 

Ug(s)+ [ ^*(U;(u))&u = 6. (4.22) 
Jo 

Hence we have ut < sup s<T Ug (s) < 00. □ 
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Proof of Theorem^ Recall from earlier remarks that it suffices to prove (14.161) . Putting 
Lemma[3]and Lemma H] together it thus suffices to show that when gt- s {') — <K w /( - > t — 
sj), with 0(A) = — ip*'(0), we have that exp{— cu(x, t)} is the unique solution to 

(I3.13p . For this to be the case, it is enough that 

-A^e-^^) - \*<f>{u}{; t - s^e-^^ 

= ^*(-\* e -<' t -^ + u* f (; t-s))- ip*(u* f {-,t - s)). (4.23) 

Note that in order to appeal to Lemma H] above, we require that g s (y) is bouned on 
each finite time horizon of s. This follows by virtue of the fact that is a Bernstein 
function (and therefore concave) and that, as indicated in the proof of Lemma HI for 
each fixed T > 0, < sup s<T sup^^d u*Ay, s) < oo. To prove (14.231) . note that 

H t _ s {-, -A*e-^~ s )) - \*<f>{u}{; t - s))e-"^ 

= X ^(,t- S )(-A*e-^-^) - \*e-»^[r'(u}(; t-s))- ^*'(0)] 

and the desired equality follows by Lemma [5j □ 

Note that up until this point in our reasoning, there are only two points where 
we have used the assumption that — ^'(0+) < oo. The first place occurs at (12. ip 
where classical literature imposes the aforesaid assumption as a sufficient condition to 
guarantee that a unique non-negative solution exists. The second place occurs is in 
justifying the finiteness in ( 14.211) . See in particular Lemma |5j 

Morally speaking the imposition of — ip'(0+) < oo in these two cases boils down 
to the same issue of allowing the application of Gronwall's Lemma to establish the 
existence of a unique non-negative solution to an integral equation. Moreover, it seems 
difficult to see how one might remove this condition in general. To see why consider, 
for example ( 14 .211) . Let / be a compactly supported function and consider all y outside 
of the support of /. It is then clear by definition that for this / and all such y, 

u* f (y,0) = 0, 

and hence Xu* } (y,o)(ty = Xo(A) = ^*(A). If in addition -0'(O+) = — oo, the failure of the 
function xo to be Lipschitz on [— A*,0] prevents us from deducing ( I4.2ip . 

We conclude by reviewing the above arguments when spatial motion is disregarded 
(formally V corresponds to a particle remaining stationary at a point) thereby giving 
the proof of Theorem |3j 

Proof of Theorem^ Suppose now that — ?//(0+) = oo and V corresponds to a particle 
remaining stationary at a point. It suffices to show that the two points (noted in the 
discussion above) where the condition — ip'(0+) < oo was used no longer need this 
assumption. 

With regard to the use of ( 12. ip . recall that, in the current setting where X t = \ \X t \\ 
for all t > 0, (12. ip collapses to the integral equation (I4.22p . Moreover, as alluded to in 
the discussion at the beginning of Section 12. 1} ( I4.22p always has a unique non-negative 
solution even when —^'(0+) = oo; see also Remark [3j 
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With regard to justifying the finiteness in (14.2 note that in the current setting, 
the quantity u*f(y, s) can be replaced by Ug(s) for 9 > 0. Since u T := inf s < T [/g(s) > 
then the estimate in (14.201) can be replaced by 

sup_ sup \xu*(-u(y)) - Xu*(-v(y))\. 

u T <u*<u T y<=R d 

Now proceeding with the proof of Lemma HI taking note of the conclusion of Lemma 
El we see that the condition — ^(0+) < oo is no longer necessary. 

□ 

4.1 Proof of Lemma [5] 

It is a strightforward algabraic exercise to deduce that 
^\X + u)-i)*{u) 
= \(^2(3u + a* + ^ (1 - e-^>e- A * x n(dx)^ 

+/3A 2 + f (e~ Xx - 1 + Xxl {x<1} )e^ x * +u ^ x U(dx). 

J(0,oo) 

It follows, with the help of IQI) . that 

= x«(A) - A^'(A*) - A / xe- {x * +u)x Yl(dx) 

</[l,oo) 

+A (2/3u + a* + J (1 - e- ux )xe- x * x U(dx)^j 
= Xu(X) + A Ufa + a*+ [ (1 {X<1} - e- ux )xe- x ' x U(dx) - ^(0)] 

V ./(O.oo) / 

However, it is again a simple exercise to deduce from (12.4p that 

ip*'(u) = 2f3u + a* + [ (l {x <i } - e- ux )xe~ x * x U(dx) 
J (0,00) 

and hence the first part of the lemma follows. 

Next notice that x u (A) is the Laplace exponent of a spectrally positive Levy process 
and therefore strictly convex and infinitely smooth on (—A*, 00). Moreover, remember- 
ing that = ip(\ + A*), we have for all A > —A*, 

X ' U (X) = ip'(X + A* + u) - ^'(A* + u) + <//(A*). 

The proof of the remaining parts of the lemma are now straightforward. □ 
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